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Abstract 

At impact energies ~ 1 GeV/u the projectile-electron excitation and loss occurring in collisions 
between highly charged ions and neutral atoms is already strongly influenced by the presence of 
atomic electrons. In order to treat these processes in collisions with heavy atoms we generalize 
the symmetric eikonal model, used earlier for considerations of electron transitions in ion-atom 
collisions within the scope of a three-body Coulomb problem. We show that at asymptotically 
high collision energies this model leads to an exact transition amplitude and is very well suited to 
describe the projectile-electron excitation and loss at energies above a few GeV/u. In particular, 
by considering a number of examples we demonstrate advantages of this model over the first Born 
approximation at impact energies ~ 1-30 GeV/u, which are of special interest for atomic physics 
experiments at the future GSI facilities. 

PACS numbers: 34.10+x, 34.50.Fa 



1 



I. INTRODUCTION 



During the last two decades there has been performed a large number of experimental in- 
vestigations of projectile-electron excitation and loss in collisions between relativistic highly 
charged ions and solid and gaseous targets. In particular, a variety of very heavy projectiles 
with the net charge 52 - 91 a.u. were used in the experiments. The experiments have also 
covered a very large interval of iinpact energies ranging from comparatively low relativistic 
energies of - 100-200 MeV/u Q-Q to extreme relativistic energy of 160 GeV/u j6|-j8| where 
the projectile velocity already only fractionally differs from the speed of light c = 137 a.u.. 

Most of the data, however, have been collected for impact energies not exceeding a few 
hundreds of MeV/u. For impact energies above 1 GeV/u there exists just a few experi- 
mental results. They include the data on the electron loss from 10.8 GeV/u Au''^"''(ls) ions 



penetrating solid targets (9|-[l0| and on the electron loss from 160 GeV/u Pb^^"^(ls) ions 
colliding with sohd p and gaseous 0] targets. 

Besides, one should also note that the experimental data on the elementary cross sections 
in collisions at very high energies collected using solid targets are not very accurate. The loss 
cross sections, reported for 160 GeV/u Pb^^"'~(ls) projectiles penetrating solid and gas targets 
differ between themselves roughly by a factor of 2 (see j6|-[3]). The reason for this difference, 
as was recently explained in [ll|, lies in multiple collisions suffered by the projectiles when 
they move in solids which does not allow one an accurate experimental determination of the 
values for elementary ion-atom cross sections. In the case of 10.8 GeV/u Au'^'^"''(ls) ions the 
data of joj and [lo[, which were both collected for collisions with solid state targets, also 
differ by about of a factor of 2 (but the reason for this is not completely clear). 

Concerning difficulties in the theoretical description of the ion-atom collisions, there 
should be mentioned two main points, which complicate the treatment of the projectile- 
electron excitation and loss: The presence of atomic electrons and the fact that in case of 
collisions with a heavy atom the atomic field may be too strong making it necessary to go 
beyond the first Born approximation. 

The latter point is known to be of extreme importance when considering the projectile- 
electron excitation and loss processes occurring in collisions of very highly charged ions 
with heavy atoms at relatively low impact energies 0.1 - 0.2 GeV/u, where the difference 
between experimental data and first Born calculations reaches an order of magnitude. At 
such energies, however, the atomic electrons play only a minor role in the projectile-electron 
transitions and their presence can simply be neglected (l3| . 

With increase in the impact energy the role of the higher order effects in the projectile- 
target interaction diminishes. Nevertheless, even at energies well above 1 GeV/u the first 
Born approximation tends to substantially overestimate transition probabilities and the 
accuracy of first Born results for cross sections remains unclear, especially in the case when 
the collision causes more than one electron of the projectile to undergo transitions. Besides, 
at such impact energies even for the most highly charged projectiles (like e.g. hydrogen-like 
uranium ions) the influence of the electrons of the atom on the projectile-electron transitions 



can no longer be ignored [13 



In the present paper we shall demonstrate that the so called symmetric eikonal (SE) 
model, extended to account for the presence of atomic electrons, can be used for treating the 
projectile-electron excitation and loss in collisions with heavy atoms at impact energies ~ 1 
GeV/u. In particular, we will show that in the limit of asymptotically high impact energies 
the transition amplitude, derived in the SE model, coincides with the amplitude obtained in 
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the light cone approximation which means that the SE model provides essentially an exact 
solution for the problem of projectile-electron excitation and loss in extreme relativistic 
ion-atom collisions. 

The forthcoming upgrade of the heavy-ion facilities at GSI (Darmstadt, Germany) will 
allow one to perform extensive experimental explorations of the different aspects of heavy 
ion - atom collisions at impact energies in the range 1-30 GeV/u. It will be seen below that 
in this range of impact energies the SE model represents a very valuable tool for describing 
the projectile-electron excitation and loss in collisions with very heavy atoms. 

Atomic units are used throughout unless otherwise stated. 



II. THEORY 

Depending on whether or not the electrons of the atom are 'active' in the collision, one 



normally distinguishes two atomic modes (see e.g. 12|, 13|) which can contribute to the 
projectile-electron transitions. In one of them, which is often called screening, the electrons 
of the atom act coherently with the atomic nucleus screening (partially or fully) the field of 
the latter. The other mode, in which the atomic electrons actively participate undergoing 
transitions, is termed antiscreening. In this mode the behavior of the projectile-electron is 
influenced mainly by the electrons of the atom while the nucleus of the atom plays only a 
minor role 

In collisions with heavy atoms, which are of special interest for the present article, the 
antiscreening mode is much weaker than the screening one and below will not be considered. 
In the screening mode the field of the atom can be regarded as external that enables one to 
reduce a many-electron problem of the ion-atom collision to a problem of the motion of the 
electron of the projectile in two external fields: the field of the nucleus of the ion and the 
field of the atom. The latter is taken as a superposition of the field of the atomic nucleus and 
the field of the atomic electrons whose space distribution is assumed to be 'frozen' during 



the very short collision time [13 



In the rest frame of the ion the electron is described by the Dirac equation 

^(r 
'di 



.'-^^mr.ty (1) 



The Hamiltonian H reads 

H = Ho + W, (2) 

where 

Ho = ccfp + f3c^ -— (3) 

r 

is the Hamiltonian for the electron motion in the field of the ionic nucleus and 

# = A(r,t)-$(r,t) (4) 

is the interaction between the electron and the field of the atom. In the above expressions 
oc and /3 are the Dirac's matrices, r coordinates of the electron with respect to the ionic 
nucleus, $ and A are, respectively, the scalar and vector potentials of the electromagnetic 
field of the incident atom. 
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The field of tfie atom in its rest frame is described by the scalar potential which, using 
results of 15 and 16 , can be taken as 



ZA4>{r') 



where 



^Ajexp(-Kjr') 



(5) 



(6) 



with the screening parameters Aj {J2j Aj = 1) and kj given in [15|] and [16|. We assume 
that in the rest frame of the ion the atom moves along a classical straight-line trajectory 
R = b + vt, where b = {bx, by) is the impact parameter and v = (0, 0, v) the atomic velocity. 
Using Eqs.([5])-([6]) and the Lorentz transformation for the potentials we obtain that in the 
rest frame of the ion the potentials of the atomic field are given by 



$(r,t) 



A(r,t) 



^ Aj-exp -KjJ-f'^{z - vty + (r^ - b)' 



V 

0,0,-$ 



(7) 



where 7 is the coUisional Lorentz factor and r = (r_|_, z) with r_|_ ■ v = 0. 
Within the SE model the transition amplitude is approximated by 



a/i(b) 



dt{xf{t) I {H-id/dt)xr{t)) 



(8) 



where the initial and final states of the electron, whose motion in the field of the ionic nucleus 
is affected by the field of the atom, are chosen according to 

dt'^t') 



00 



(9) 



Here, ipo and ipn are the initial and final undistorted states of the electron in the ion. 

Making use of the fact that the dependence of the scalar potential $ on the electron 
coordinates and time is of the form $ = 'yZAf{s±,^\z — vt\), where s_|_ = r_|_ — b the 
transition amplitude can be transformed into 



a/i(b) 



c 
I— 



V J-00 

(m 

[ 1' 



dtexp{iUnot) [iJn 



exp (^i dt'^it')^ X 



(10) 



where a;„o = ^n — ^o is the electron transition frequency and V_l denotes the two-dimensional 
(in the (x, ?/)-plane) gradient operator. 

The amplitude ( JTOl) can be simplified by employing the relation 



lim 



rftexp(icj„o^) exp(-A|t|) / dt'^{t') 

J —00 



/ dtexp{iuJnot)^{t), {(^nO 7^ 0), 

UlnO -J -co 



(11) 
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that yields 



a/i(b) 



I— I dt exp{iUnot) (i^n 

V J-oo 



V 



a. 



r 



I (Vx$(t)) - ax 



exp j dt'<^{t')J X 



(12) 



The Limit of weak interaction 



If the interaction between the electron of the ion and the atom is sufficiently weak one 
can replace the exponent exp (j, dt'^{t')^ in (fT2l) by 1. Then integrating in (fT2l) by parts 
over the space and using the continuity equation 



dp, 



nO 



dt 



■ JnO 







for the transition charge and current densities, 



Pno = 'ipi ^0 exp (iunot) 

jnO = COL tpo exp (iUnot) 



we obtain 



«/i(b) 



i- dt expUuJnot) ( Ipn 

V J-oo \ 



V (9$ 



^0 



Taking into account that ^ = and 



(9$ r+oo 

dt exp (iunot) = -iuJno / dt exp (iunot) $ 

ut J —oo 



we arrive at the transition amplitude 

r+co 

a/i(b) = i / (itexp(za;no^) (T/'n 



$ ( 1 - -a, 



(13) 



(14) 



(15) 



(16) 



(17) 



which coincides with expression for the amplitude obtained in the first Born approximation. 



B. The High-Energy Limit 

An important question, which will be addressed in this subsection, concerns the high- 
energy limit (7 — 00) of the SE model. Keeping in mind that $ = ■-fZAf{s±,'j\z — vt\) one 
can show that 

dtexp{iUnot)^{t) = exp{iUnoz/v) / d^f{s±,^) 



exp {iunoz/ v) G isj_, — I , (18) 
V 1^ J 
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where at the moment the exphcit form of the function G is not important. Correspondingly, 



/ + 00 
dt =G(sx,0) =Gq{sa_). 
-oo 



(19) 



At sufficiently large impact energies, where the difference between G = G (^s±, and Go 
essentially vanishes, we can replace the amplitude ( fT2l) by the following expression 



af^(h) 



. c 



(V'n |exp (iG) exp{iuJnoz/v)Ga;, \ tpo) 



- I {ilJn |exp {iG) exp(iw„o^/f) (V^G) ■ Q:±| ^o) 



(20) 



We now take the second line of fl20|) . integrate there by parts, use the continuity equation 
f lT3|l and then again integrate by parts. As a result of these manipulations, expression f l20|) 
transforms into 



a/i(b) 



exp(za;„o^/^) exp (zG) ^1 



-a, 



4^0 



{ipn \exp{iujnoz/v) exp (iG) (1 - iG) V'o) 



The limit 7 — ?■ oo of the amplitude (1211) is given by 



afi{h) 



ey.Y>{iu]noz / v) exp {iGq) 



-a, 



(21) 



(22) 



and it coincides with the transition amplitude derived in the so called light-cone approach 



see 



The light-cone approach is strictly valid at 7 — ?■ 00 and in this limit enables one to solve 
the problem of electron loss (ionization) and excitation exactly. One should also mention 
that the case of a strong ion-atom interaction this exact solution does not coincide with 
the first Born results, no matter how high is the impact energy [l3|. Taking the above two 
points into account we can make the following conclusions. First, in the high-energy limit 
the SE model yields an exact solution for the transition amplitude of the projectile-electron 
excitation and loss. Second, even at 7 — 00 the results of the SE model still in general differ 
from those of the first Born approximation. 

It is rather obvious that such conclusions would also hold if the electron transitions in 
the ion are caused by the collision with a charged particle (e.g. a stripped atomic nucleus) 
and, for instance, can be applied to /T-shell ionization of atoms by high-energy bare nuclei. 
In this respect one should note that in the literature on relativistic ion-atom collisions there 
had been already attempts to consider the high-energy limit of distorted-wave models for 
the case of atomic ionization or excitation by the impact of a nucleus, where these processes 
are treated as a three-body Coulomb problem (the incident and atomic nuclei and atomic 



electron). In particular, starting with the work of [18|-[20|, it had been assumed that the 



high-energy limit of such distorted-wave models for ionization and excitation processes, like 

and the 
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the continuum-distorted-wave-eikonal-initial-state (CDW-EIS) (see e.g. 
SE, is simply that of the first Born approximation. 

We have just seen, however, that for the SE model this is not true and that in the high- 
energy limit the transition amplitude, obtained in this model, goes over into the amplitude 
derived in the light-cone approach. Since at 7 3> 1 the CDW-EIS model becomes essentially 
identical to the SE, the high-energy limit of the CDW-EIS coincides with the light-cone 
approach but differs from the first Born approximation. 
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C. The explicit form of the amplitude 



Using the expUcit form ([7]) of the scalar potential we obtain that 
dtexp{iUnot)^{t) 



exp (iujrioz/v) J2 ^j^o {s±Aj) 



V 



where where Kq is the modified Bessel function [23j, s_l = |sx| = |rx — b| and 
Besides, it also follows from (l23l) that 



2Z, 



Taking Eqs. (l23l) and (l25l) into account the transition amplitude becomes 



.2Zac ^ I 



exp \ 2 "° ^ exp [ i — — ^ A^Kq [k^ I x 



V J 



^""Kq {s±Aj) - i ^'^^ Ki (s_lAj) s_l ■ q;j 



.7 



4^0 



where Ki is the modified Bessel function 23 . 



(23) 



(24) 



(25) 



(26) 



D. The limit of vanishing screening 

The results obtained above, can also be applied for treating ionization and excitation of 
neutral atoms in relativistic collisions with bare nuclei. This can be done by setting in fl26|) 



Kj = 0, replacing there by Z^r, where Z^ is the charge of the nucleus incident on the 
atom, and regarding now ipo and ipn as the initial and final states of the 'active' electron of 
the atom. Keeping in mind that J2j Aj = 1, taking into account that Ko{x) ^ — In (^|^ — F 
for I X 1^ 1 (see e.g. [13]), where F is Euler's constant, and disregarding an inessential 
phase-factor we obtain 



a/i(b) 



.2Znc I , 



UnOZ 

exp I I I exp 



.2Za, 
-I m s_L 

V 



— Kq S± - I 



Ki (s_lA,-) s_l ■ q;_l 



(27) 



III. APPLICATIONS 

The SE model is normally regarded as a tool for describing collision-induced transitions 
between bound states. In our case it would mean that this model should be first of all 
applied for treating the projectile-electron excitation. It is known that at compara- 
tively low relativistic impact energies the model has a problem with describing transitions 
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FIG. 1: Cross section for the electron loss from Au^'^''~(ls) ions in collisions with neutral Au atoms. 
Solid curve: eikonal results. Dash curve: first Born results. The circle displays experimental data 
from 0] while the square shows the result of 10 1 scaled to the gold target. 



involving electron spin flip. However, this problem diminishes when the energy increases 
and practically disappears at energies ~ 1 GeV/u. 

The SE model can also be used for considering the projectile-electron loss (in particular, 
the total cross section). In such a case the model also becomes more accurate when the 
impact energy increases. In particular, provided in the rest frame of the atom the magnitude 
of the velocity of the electron, emitted by the projectile, is of the order of the collision 
velocity, the SE model can be applied for calculating not only the total but also differential 
loss cross sections. In practical terms this condition holds starting already with 7 ~ 2-3, 
i.e. at impact energies ~ 1-2 GeV/u. 

As was shown in the previous section, at asymptotically high impact energies the SE 
model yields an exact solution. In the case of the projectile-electron excitation and loss in 
collisions with neutral atoms even for the most highly charged projectiles the region of such 
energies is actually reached already at 7 ~ 30-50. Thus, the model should work excellently 
starting with the magnitude of 7 of a few tens. 

Taking all this into account one can expect that the SE model performs quite well at 
impact energies ~ 1-30 GeV/u, which are of special interest for the present study. Below 
the projectile-electron excitation and loss will be considered in this energy range for col- 
lisions between hydrogen- and helium-like highly charged ions and heavy atoms by using 
the SE model. Results of this model will also be compared with those of the first Born 
approximation. 



A. Single-electron loss 

In figure [U we show results for the electron loss from incident Au^^"''(ls) projectiles in 
collisions with neutral Au atoms at impact energies 1-30 GeV/u. These results include our 
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first Born and eikonal calculations as well as experimental data from [9| and 10 1 on the 
electron loss from 10.8 GeV/u Au''^^(ls) projectiles. 

In the case considered our eikonal and first order results differ by 15-35%. As expected, 
when the impact energy increases the difference between them decreases. On overall, the 
difference is not large but nevertheless should be taken into account if precise cross section 
values are needed. 

Both the first Born and eikonal cross sections agree neither with the experimental data of 
joj nor with that of [lo[. The data from joj are substantially smaller (by about of 30%-50% 
) while the data reported in 10| are considerably larger (by about a factor of 1.3-1.4) than 
our results (25j . 

As was already mentioned in the Introduction, there is a difference by roughly a factor 
of 2 between experimental cross sections reported for for the loss from 160 GeV/u Pb^^"*" 
ions in collisions with solid and gas targets with the solid state cross sections being larger. 
The origin of this difference, which was explained in [ll|, lies in multiple collisions between 
the projectile and atoms inside solids which effectively enhance the electron loss process. 
Since the magnitude of this difference depends on the impact energy and decreases when 
the energy decreases, similar reasons are probably responsible for the observed disagreement 



between our results for atomic targets and the experimental data from [10 



B. Two-electron transitions 



If a projectile- ion initially carries several electrons, then more than one electron of the 
ion can be simultaneously excited and/or lost in a single collision with a neutral atom. 
Helium-like ions are simplest projectiles, for which simultaneous transitions of more than 
one electron are possible, and below we consider double electron loss and loss-excitation for 
the case of such projectiles. 

It is known (see e.g. [l3j) that, provided the condition '^'^^ > 0.4 is fulfilled, two- 
electron transitions in a heavy helium-like ion occurring in collisions with an atom are 
governed practically solely by the independent interactions between the atom and each of 
the electrons of the ion. In order to describe such transitions one can apply the independent 
electron model. According to this model the cross section for the double electron loss from 
a helium-like ion is given by 

ai,i = 271 / dbbPi,iib), (28) 
Jo 

where the probability Pij{b) for the two-electron loss is given by 

PiAb) = PLsib), (29) 

where Piossib) is the single-electron loss probability in a collision with a given value of the 
impact parameter b. 

The cross section for the simultaneous loss-excitation in the case of a helium-like ion is 
evaluated as 

ae,l = 271 / dbbPe,i{i>). (30) 
where the probability P{b) for the two-electron process is given by 

Pe,m='^Pexc{b) Piossib), (31) 

where Pexc{b) is the single-electron excitation probability. 
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FIG. 2: Cross sections for double electron loss from Au^'^"*"(ls^) ions in collisions with neutral Au 
atoms. Solid curve: eikonal results. Dash curve: first Born results. 

1. Double- electron loss 

Figure [2] shows calculated cross sections for double electron loss from Au^^'^(ls^) projec- 
tiles incident on neutral Au atoms at impact energies 1-30 GeV/u. Compared to the single 
loss, now the difference between the eikonal and first Born results is much more pronounced. 
For impact energies ~ 1 - 5 GeV/u these calculations predict even qualitatively different 
dependencies of the cross sections on the collisions energy. The absolute difference between 
the first Born and eikonal cross sections ranges between ~ 2.5 at 1 GeV/u to ~ 1.5 at 30 
GeV/u. Moreover, as additional calculations show, for impact energies above 30 GeV/u the 
ratio of ^ 1.5 remains almost a constant and, thus, even at asymptotically high collision 
energies the first Born calculation still substantially overestimates the cross section values. 

2. Simultaneous loss- excitation 

In figure [3] we present results of calculations for the simultaneous projectile-electron ex- 
citation and loss, U^°+(ls^) — U^^+(n = 2,j) + e~, where n and j are, respectively, the 
principal quantum number and the total angular momentum of the excited state of the 
hydrogen-like ion. The projectile is assumed to collide with Kr, Xe and Au atomic targets 
at an impact energy of 20 GeV/u. Similarly to the case of the double electron loss, con- 
sidered above, we observe that the differences between the cross sections, calculated with 
the first Born and eikonal probabilities, can be quite substantial if the atom is sufficiently 
heavy. For collisions leading to the population of the states with j = 1/2 the first Born 
results are larger by a factor of ~ 1.16 (Kr), ~ 1.38 (Xe) and ~ 1.84 (Au). When the states 
with j = 3/2 are populated this ratio is ^ 1.12, ~ 1.31 and ^ 1.67, respectively. 

We see that the difference between the first Born and eikonal cross sections turns out to be 
somewhat smaller for transitions involving the states with j = 3/2. This can be attributed 
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FIG. 3: Cross sections for the simultaneous electron loss-excitation, U^'^"'"(ls^) — )■ U^^~^(n = 2,j) 
+ e~, occurring in collisions with Kr, Xe and Au atomic targets at 20 GeV/u. j = 1/2 and j = 3/2 
are the angular momentum of the states of the hydrogen-like uranium ion. Squares and circles 
show results for j = 1/2 and j = 3/2, respectively. Open and solid symbols denote the cross 
sections obtained using the first Born approximation and the SE model, respectively. 
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FIG. 4: Cross sections for excitation of 20 GeV/u U^^"*"(ls) projectiles into the states with n = 
2,j = 1/2 (squares) and n = 2,j = 3/2 (circles) in collisions with Kr, Xe and Au atomic targets. 
The cross sections were calculated using the SE model. 



to the fact that, compared to the j = 1/2 case, these transitions are characterized by larger 
impact parameters. As a result, the field of the atom acting on the electrons of the ion is 
weaker and, therefore, the first Born treatment becomes less inaccurate. 

As additional calculations show, the differences between the results of the eikonal and first 
Born calculation does not substantially change when the impact energy increases further. 
Thus, like for the double electron loss, even at asymptotically high impact energies the 
first Born calculation may considerably overestimate the cross section for the simultaneous 
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FIG. 5: Probabilities for the electron loss, excitation to the states with (n = 2, j = 1/2) and 
(n = 2, j = 3/2) in collisions of 20 GeV/u U^^"^(ls) ions with Au atoms. 

projectile-electron excitation and loss. 

According to both the first Born and eikonal results the cross sections for the simultaneous 
loss-excitation are larger for transitions to the states with j = 1/2. However, the pure 
excitation U9°+(ls) ^ U^°+(ls; n=2, j) (or the excitation U9i+(ls) -> \J^^+{n=2, j), see 
figure H]) at these rather high energies proceeds already more efficiently into the states with 
j = 3/2. The origin of this interesting peculiarity can be understood by considering the 
single-electron transition probabilities. It turns out that the probability for the electron loss 
has a larger overlap with the probability for the excitation into the states with j = 1/2 (see 
for illustration figure E]) . 

IV. CONCLUSIONS 

We have considered the symmetric eikonal model for treating projectile-electron transi- 
tions in relativistic collisions with heavy atoms. 

We have shown that at asymptotically high impact energies this model yields the same 
results as the light-cone approach and, thus, offers an exact solution for the transition ampli- 
tude. In the limit 7^1 another popular distorted-wave model - the continuum-distorted- 
wave-eikonal-initial-state approximation - becomes essentially identical to the symmetric 
eikonal model. Therefore, at asymptotically high impact energies the CDW-EIS also be- 
comes exact and, contrary to what was stated earlier, does not coincide with the first Born 
approximation. 

Using the SE model and the first Born approximation we have calculated cross sections 
for the projectile-electron excitation and loss. We have shown that at impact energies 
~ 1-30 GeV/u, which are relevant for the future GSI facility, there are noticeable deviations 
between the eikonal and first Born results for single-electron transitions and very substantial 
differences between such results for transitions involving two projectile electrons. Moreover, 
these very substantial differences 'survive' even in the asymptotic limit 7 — )■ 00. 

Our results clearly show a great advantage of the SE model over the first Born approxi- 
mation. 
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